Fourth— Order Ricci Gravity 



A. Borowiec a , M. Francaviglia 6 and V.I. Smirichinski c 

a Institute of Theoretical Physics, Wroclaw University, Poland 

b Departimento di Matematica, Unversita di Torino, Italy 

c Joint Institute for Nulclear Research, Dubna, Moscow Region, Russia. 

Abstract 

The Euler-Lagrange equations of motion for the most general Ricci type gravita- 
tional Lagrangians are derived by means of a purely metric formalism. 

Einstein metrics are extremals of the Einstein-Hilbert purely metric variational problem. 
It is known that the non-linear Einstein-Hilbert type Lagrangians f{R)y^g, where / is a 
smooth function of one real variable and R is the scalar curvature of a metric g [], lead 
either to fourth order equations for g which are not equivalent to Einstein equations unless 
f(R) = R — c (linear case) or to the appearance of additional matter fields. 

Fourth-derivative theories of gravity enjoy better renormalizability properties then stan- 
dard gravity based on the Hilbert-Einstein Lagrangian density. Besides, higher-order the- 
ories are of interest in cosmology, since they do contain every vacuum solution of GR 
(including Schwarzschild solution). Linear Ricci type Lagrangians have been studied in 
this context by many authors (see e.g. [1, 7-10] and references therein). 

In the present note we are going to use variational decomposition techniques (see e.g. 
[4]) in order to obtain equations of motion for the most general (nonlinear) Lagrangian 
densities of the Ricci type (see formula (6) below) by means of a purely metric formalism. 
As an intermediate step one performs the variational decomposition in the framework of 
the so-called Palatini formalism, i.e. keeping a metric and a (symmetric) connection as 
independent variables [2-6]. 

In the sequel we shall use lower case letters and rp v = r^ au to denote the Riemann 
and Ricci tensor of an arbitrary (symmetric) connection T: 

a _ a /p\ _ rs -pa _ □ p« , pa per _ pa pff 

i.e. without assuming that V is the Levi-Civita connection of a metric g. We shall assume 
that space-time is four- dimensional but our results will be valid in any dimension n > 3. 

As it is also known, the linear "first order" Lagrangian r^/g, where r = r(g,T) = 
g af3 r a p(T) is a scalar concomitant of a metric g and a linear (symmetric) connection F, f\ 
leads to separate equations for g and F which turn out to be equivalent to the Einstein 

1 One simply writes yfg for y^detg]. 

2 Now, the scalar r(g,T) = g a ^r a p(T) is no longer the scalar curvature, since T is no longer the 
Levi-Civita connection of g. 



equations for g (so-called Palatini principle, c.f. [1-4]). In other words, the first-order 
metric-affine Lagrangian r^fg and the purely metric Einstein-Hilbert Lagrangian R^fg are 
equivalent in the sense that they lead to the same class of solutions, namely to vacuum 
Einstein metrics. 

Unlike in the purely metric case, an equivalence with General Relativity also holds for 
non- linear gravitational Lagrangians 

L f (g,r) = y/gf(r) (2) 

(parametrized by the real function / of one variable) when they are considered within 
the first-order Palatini formalism. More exactly, it was shown [6] that (with the notable 
exception of some special non-generic choice of the function /) we always obtain Einstein 
equations as gravitational field equations. This property of Einstein equations was called 
universality. 

The universality property holds also true for "Ricci squared" non-linear Lagrangians 

L f (g,r) = y/gf(8) (3) 

where s = s(g, T) = g^g^r^^r^ and = r^(T) is the symmetric part of the Ricci 
tensor of T (see [3]). (Thereafter round brackets around indices denote symmetrization) . 
But in this case, besides the Einstein equations, one obtains algebraic constraints which 
can be interpreted as an almost-product or an almost-complex structure on the space-time 
manifold [2]. Moreover, which one of these two can be realised depends on the sign of a 
cosmological constant which is introduced within the calculations (see [3] for details). 

Let us consider a (1,1) tensor valued concomitant of a metric g and a linear torsionless 
connection T defined by 

S£ = S£(3,r) = ^r (JU ,)(r) (4) 
One can use it to define a family of scalar concomitants of the Ricci type 

s k = trS h (5) 

for k — 1, . . . , 4. We can in fact eliminate the higher order Ricci scalars s k with k > 4 
by using the characteristic polynomial equation for the 4x4 matrix S (see [5]). One 
immediately recognizes that r = Si — trS and s = s 2 = trS 2 . 

The most general family of non-linear Ricci-type gravitational Lagrangians in four- 
dimensions, i.e.: 

Lp(g, r) = y /gF(s 1 , s 2 , s 3 , s 4 ) (6) 

is thence parametrized by a smooth real-valued function F of 4-variables. This includes 
the Lagrangians (2) and (3) as particular cases. 

As already mentioned above we start with the "Palatini variational principle", i.e. we 
choose a metric g and a symmetric connection T as independent dynamical variables. 
Variation of Lp gives then 

5L F = ^gS g F- X -Fg afi Sg * + ^g5 r F (7) 



where obviously 5F = YX=i $ s k an d ^k = fr - We see at once that 

S g s k = fetr^*" 1 ^) = k(S k - 1 y a r(M Sg a/3 
which is clear from 5s k = ktr(S k ~ 1 8S). Accordingly one has 

5 g F = ^r (M 6g^ (8) 
where for simplicity we have introduced a (1,1) tensor field concomitant 

4 

k=l 

In a similar way one calculates 

S r F = T a a 5r (a/3) = Sr (a(3) (9) 

where, as usual, the metric g has been used to rise and lower indices. Moreover, as it can 
be easily checked by a direct calculation, the (0,2) tensor field 

-pa/3 — pa cr/3 

is already symmetric (see [5]). 

Substituting (8) and (9) into equation (7) gives 

SL F = y/g \Tl r {M - ^Fg aP ] 5g a ? - yfgF* 8r {a(3) (10) 

Now, taking into account the so-called Palatini formula 

Sr {a/3 ) = V^r^ - V {a 5T a p)a 

V Q being the covariant derivative with respect to T and performing the "covariant" Leibniz 
rule one gets the variational decomposition formula under the form 

SLf = v^[^^)-^^]^-V w [v^(^63:-^^)]^ + div (11) 

which splits 5Lp into the Euler-Lagrange part and a boundary term generically denoted 
by div since we shall not exploit its explicit form here. The boundary term, however, 
plays an essential role during the calculation of a canonical energy momentum complex 
(see e.g. [4]). This will be therefore investigated in more details elsewhere. 

Accordingly, the Euler-Lagrange field equations ensuing from the Palatini formalism 
read as follows 

^ a r {p)a) = l -Fg aP (12) 
V, [y/gF *] = 0. (13) 



Their analysis have been already performed in [5]. Also for this case the universality of 
Einstein equations takes place. 

Our goal in the present note is to calculate the Euler-Lagrange equations in the purely 
metric framework, i.e. assuming to be the Levi-Civita connection of g. For this 
purpose one has to continue our calculations with ST^ replaced in (11) by 

Waft = ^(V Q % + VpS9cuT ~ V,^)- 

For example, transvecting the first term above with the second part of (11) one calculates: 

V„ (y/g SI - 5{\) \g x °VJgp a = div - l -V a V u (y/g 5» x - 5g Pa 

Similar calculations for the two remaining terms give rise finally to the fourth-order Euler- 
Lagrange equations for the Lagrangian densities (6) which, after some calculations, can 
be presented in the following form: 

-DF*? + -# Q/3 V ' a V p T ap - V a V {a ^ )a - -g a(3 F + T a(pi R^ = (14) 
2 2 2 

where □ denotes the d'Alambertian (or Laplace-Beltrami) second order differential oper- 
ator and R^ v is the Ricci tensor of g. Since JF Q/3 depends on the second-order derivatives 
of g equation (14) provides a fourth-order non-linear operator acting on g. 

Example 1. Take F = f(R). Then T afi = f'(R)g al3 and (14) yields the well-known 
result 

[ g ^n _ v a V?]f(R) - \f{R)g^ + f'(R)K* = (15) 

Example 2. Take F = f(R^R^); then T a<i = 2f(R^R^) R a/3 . In this case (14) gives 

□ (/^)+^V CT V p (/^ CT 0-2V CT V^ = 0, (16) 

which for the linear case f — 1 reduces to 

DRaf 3 + g ap v pR a P _ 2 V a V {a (R^) - \g ap R^R^ + 2R^ a 'It® = 0. (17) 
This can be recasted (using Bianchi identities) as follows (compare [9]): 

□i? Q/3 + \g^UR - V a V p R - \g af5 R^R^ + 2R f3<Tap R ap = 0. (18) 
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